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ABSTRACT 

Searches for planetary transits carried out in open and globular clusters have yielded to date only a handful of weak, unconfirmed 
candidates. These results have been interpreted either as being insignificant, or as evidence that the cluster chemical or dynamical 
environment inhibits the planetary formation or survival. Most campaigns were limited by small sample statistics or systematics from 
ground-based photometry. In this work we performed a search for transiting planets and variables in a deep stellar field of NGC 
6397 imaged by HST-ACS for 126 orbits. We analyzed 5,078 light curves, including a pure sample of 2,215 cluster-member M0-M9 
dwarfs. The light curves have been corrected for systematic trends and inspected with several tools. No high-significance planetary 
candidate is detected. We compared this null detection with the most recent results from Kepler, showing that no conclusive evidence 
of lower planet incidence can be drawn. However, a very small photometric jitter is measured for early-M cluster members (< 2 mmag 
on 98% of them), which may be worth targeting in the near future with more optimized campaigns. Twelve variable stars are reported 
for the first time. 

Key words, techniques: photometric - stars: planetary systems - clusters: individual: NGC 6397 
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1. Introduction 

More than seven-hundred extrasolar planets are known (exo- 
planet . eu database). Most of them are characterized only by 
means of radial velocities (RV) and lack any information about 
their "real" mass Mp, because it can be measured only Mp sin /, 
being / the inclination of the orbit with respect to the line of 
sight. Their size is also unknown, which prevents us from get- 
ting any clues about their density and physical composition. 
Exoplanetary transits are highly complementary to RV tech- 
niques, providing the planetar y radius fro m the stellar radius 
Ri, in a direct geometrical way (ISeageij201 iL p. 55). Photometric 
searches for transits can also go much deeper than RVs in mag- 
nitude, and can monitor thousands of stars simultaneously. The 
Kepler mission ( Borucki et al.i 120101) . for instance, is demon- 
strating the power of this technique by discovering many plan- 
etary systems with unexpected properties. Hundreds of Kepler 
"candidate planets", for which confirmation and mass measure- 
ment via RVs re mains infeasible, are still very useful for sta- 
tistica l purposes (iHoward et alJ 120111 : ISchlaufman & LaughlinI 
I20T1I) 

Star clusters, and in particular globular clusters (GC), of- 
fer a unique opportunity to study how the chemical and dy- 
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namical environment affects the planetary formation and evo- 
lution. They are also comprised of stars that share (in most 
cases) the same age and chemistry, and whose radii Ri, 
and masses are reliably known on their main sequence 
(MS). Open clus ters (OC) have been t argeted for extensive 
transit searches (Mocheiska et alJ 120051: iMontalto et all l2007l 
2011; Hartman et al. 2009, most notably) but only a handful 
of weak, unconfirmed candidates have bee n detected so far 
(iMocheiska et al.ll2006l:lMontalto eTani201 l l). A global reanaly- 
sis suggests that the overall statistical significance of these cam- 
paigns is so low that it could be compatible with the planet host 
incidence observed in the field (van Saders & Gaudi 201 1). GCs, 
on the other hand, are on average much richer in stars than OCs, 
providing more statistical significance in case of null detection. 
The only GCs monitored for t ransits have been 47 Tucanae for 
eight days with HST WFPC2 (Gilliland et al. 2000) and for 33 
days from the ground ( Weldrake et al. 2005); and to Centaur i 
with a 25-day ground-based campaign dWeldrake et alJ 2008h. 
No pl anetary detection has been claimed, and iGilliland et al] 
(l2000h concluded that the planet occuiTence in 47 Tuc is smaller 
by a factor of ten compared to field stars. The reasons which have 
been mostly hypothesized to explain the lack of giant planets in 
GCs are their metallicity, and their dynamical environment. 

It has long been known that metallicity is a strong primary 
parameter that coii'elates with the fraction of sta rs with plan- 
ets Op. For giant planets iFischer & Valentil (l2005l) . among oth- 
ers, showed an increase from the typical value ~ 0.03 for 
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Stars with solar metallicity up to <i>p ~ 0.25 for very metal-rich 
([Fe/H] > +0.3) stars. It is still disputed whether for moderately 
metal-poor stars (-0.5 < [F e/H] < 0) this correla t ion becomes 
flat o n values Op ^ dUdrv & SantosI 120071: ISantos et al] 
l2011h . or Op conti nues to decreas e exponentially towards lower 
metalUcities (Joh nson et al.ll20l"ol) . As more and more low-mass 
planets (Mp < 30M^) are being discovered, it becomes clear that 
their Op is much larger than that of giant planets (Lovis et al, 
l2009tlWittenmver et al.ll201 lb . probably around -20-30 %. The 
occurrence of low-m ass planets seems to be insensitive to the 
host star metal licity ( Udrv et aL| 20061), except maybe for very 
late-ty pe stars dJohnson & A pps 20091: ISchlaufman & LaughUnI 
l201lh . but M-dwarf metallicity calibration is still too uncertain 
to draw conclusions. 

The complex dynamical environment of a cluster is the sec- 
ond major concern about the survival of protoplanetary sys- 
tems. Theoretical studies give different answers, but some of 
them imply that gravitational stripping processes are not enough 
strong to disrupt very short-period (P < 5 days) planetary sys- 
tems, even in the densest regions of a typical GC jFregeau et alj 
[2OO6; Spurzem et al. 2009). Indeed, planets do exist in clusters. 
Giant planets have been de tected around evolved stars be long- 
ing to open clusters (L ovis & Mavoij2007l : Sato et alj2007h. an d 
around a pulsar in the globular cluster M4 (iThorsett et al ]fl999l) . 
It is therefore difficult to understand why no planets were de- 
tected in 47 Tuc. The search by Gilliland et al. (2000) was sensi- 
tive only to giant planets (Rp > IRjup) around stars in the upper 
MS, excluding late-K and M dwarfs. More data are needed to 
sample other regions of the parameter space. 

On the other hand, it is convenient to search for transits 
around KM main sequence stars because of their larger expected 
signal. These targets, even in the nearest GCs, are faint (V > 18) 
and extremely crowded. Space-based observations are necessary 
to achieve the needed signal-to-noise ratio (S/N) per time unit, 
and to minimize the number of false positives due to blended 
photometric contaminants. The wide-field imagers mounted on 
HST (ACS and WFC3) are unrivalled in this regard, and their 
large archive of deep photometric series of GCs is already avail- 
able to exploit for transit searches. However, we will show in 
Sec.[3]that HST time-resolved photometry is influenced by many 
sources of systematic eiTors that require a careful coiTection. We 
developed specific tools to this purpose, and applied them to a 
test case. Other data sets (such as those on 47 Tuc) will be in- 
vestigated in the near future, to test whether a more optimized 
campaign is worth pursuing. 

In this Paper we present a search for variables and plan- 
etary transits in NGC 6397, exploiting a 126-orbit data set 
from the HST Advanced Camera for Surveys (ACS: GO- 
10424; PLRicher). NGC 6397 is a metal-poor ([Fe/H] ^ -2) 
core-collapsed globular cluster, and the second-closest known 
dGratton etal J 120031; [Richer et al]|2008h . Our work is in some 
sense complementary to that presented by Stello & Gillilai^ 
(2009), who employed the same data set to study the microvari- 
ability of highly saturated red giants. Though these data were not 
optimized for a transit search, its unprecedented depth allowed 
us to search for planets on a homogeneous sample of 2,215 
member M-dwarf stars down to the hydrogen-burning limit. M 
dwarfs are considered the most promising targets to discover 
rocky planets in the near future tScalo et alJl2007,) . They are the 
smallest stars (7?* ^ 0.1-0.5/?©), and therefore the transit depth 
for a given planet is increased by a factor of ~ 4-100 com- 
pared to a solar-type host star They are also intrisically faint 

^ 0.02-5 ■ IO^^Lq), meaning that an habitable planet would 
have an orbital period of only P = 10-30 days. 



2. Observations and data reduction 

Our analysis is based on the HST data set GO- 10424 (PI: 
Richer), originally aimed to probe the bottom of the MS and 
the end of the white dwa rf (WD) cooling sequence in NGC 6397 
(Ha nsen et all 2007; Ri cher et al.ll2008h . A single 202" x 202" 
field, located 5' from the cluster core was imaged for 126 orbits 
with the wide-field channel of the ACS. In each orbit a single ex- 
posure through the F606W filter was bracketed by two F814W 
images, with exposure times ranging between 584 and 804 s 
(median: 704 s). Overall, 252 F814W frames and 126 F606W 
frames were secured, for a total exposure time of ~ 50 hours. 
Shorter exposures (fexp - 1, 5, 40 s) were also taken to measure 
the brighest stars, but they have not been used in our study. 

The dynamic range of the "deep" exposures is perfect to our 
purposes, as the observed luminosity function (LF) of the cluster 
MS members peaks at wpguw - 21 (that is, on well-measured 
stars with S/N ~ 200 in the same filter; Fig. [T] histogram on 
the right panel). Saturation occurs around MOV spectral type 
at mpsMw - 18.7 and the faintest cluster members (M9V) are 
found at mf^i4w - 24 (S/N ~ 20), so that a sample of ~ 2,000 
cluster M dwarfs down to the hydrogen-burning limits is avail- 
able to transit search (Fig.[T] left panel). 

On the other hand, these data were not optimized to search 
for transit events, therefore two aspects of the observing setup 
are somewhat limiting. Fi rst, the time coverag e of the frames is 
not continuous, as it is in iGiUiland et al.l (j2000). The ~ 50 h of 
integration time were split among twenty-one non-contiguous 
days, spanning a twenty-eight day period. This translated into 
a much lower completeness for our search, especially for the 
longest periods (P > 3 d, see Sec.|5]l. The second reason is that 
each pointing was shifted with a ten-position dithering pattern 
Ax, Ay plus a subpixel offset 6x, 6y. While u sually dithering is 
a winning choice on undersampled images dAnderson & Kind 
2000), it prevented us from reaching the best photometric ac- 
curacy possible on the brightest stars. On these targets, the 
amount of random noise is so low that the unavoidable flat-field 
and pixel-to-pixel residual errors become no more negligible. A 
zero-point correction, discussed in next Section, was developed 
to suppress these systematic errors. 

We carried out the data reduction on individual bias- 
corrected and flat-fielded . fit images provid ed by the HST 
pipeline. We employed the master input list from Anderson et alj 
(2008), and the code described by Anderson & King (200^ 
bas ed upon the effec tive PSF (ePSF) approach first developed 
bv Anderson & Kinj (12000) . Four tests have been performed on 
a subset of twenty F814W frames taken at the same integer-pixel 
dithering position, to choose the best reduction strategy between: 

1 . allowing spatially-consta nt perturbed ePSFs as described by 
lAnderson & KingI (l2006h : 

2. correcting the raw frames for charge transfer (in)efficiency 
(CTE) w ith the pixel-based algorithm proposed by 
Anderson & BedinI ( i2010.) ; 

3. using both the 1) and 2) corrections; 

4. using neither. 

For all the measured sources, the RMSs of their light curve as a 
function of magnitude were compared. Methods 1-3 above pro- 
vided no detectable improvement over the fourth choice, there- 
fore we decided to not apply any CTE or PSF correction at this 
stage. In fact, actual PSFs are both spatially and temporally vari- 
able, and require a frame-to-frame a posteriori correction that 
will be explained in Section[3] 
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instrumental magnitude F814W 



Fig. 1. Left panel. Color-magnitude diagram (mpeoew - 'Wf814w, 'Wf814w) for all the stars in the 'Anderson et alJ (l2008h master list; 
the 5,078 sources selected in this wor k lie betw een the saturation limi t at wpgHw - 19 an d the faint limit ;7Jf8i4w - 26.5 (dash-dot 
lines). The red line is the isochrone bv lDotter"et al. (2007) employed in lRicher et al.l (l2008h . The dotted line marks the loci occupied 
by equal-mass MS-MS binaries. Right panel. RMS for light curves in our sample as a function of the instrumental magnitude. 
Red circles: median RMS averaged over 0.5 mag bins, without any correction. Cyan circles: the same, after PSF-ZP correction. 
Green circles: both PSF-ZP and dither-ZP corrections applied. Small points are the RMS of the individual light curves after both 
corrections. The solid orange line is the expected theoretical noise level. The superimposed histogram represents the number of 
targets in each 0.5 mag bin (scale at the right). 



Hereafter we focus our analysis only on F814W frames 
for many reasons. They are much deeper than F606W im- 
ages on faint, red, low MS stars, and are less affected by PSF 
short-term variations as showed for this very same dataset by 
[Anderson & Kind (l2006h . The sky background is also much 
lower in F814W exposures. The way the F814W frames were 
sampled (at the beginning and the end of each 96-minute orbit) 
makes F606W frames nearly useless to increase the transit de- 
tection efficiency. In this way we also avoided the need to correct 
for tricky registration of the light curves between the two filters. 

Sources beyond the saturation limit on the longest exposures 
(mpsuw ^ 19) and sources detected on less than 200 frames 
(mF8i4w ^ 26.5) have been excluded from this study, leaving 
5,078 objects which include cluster members, field stars, and a 
limited number of non-stellar objects. We evaluated the mem- 
bership of each entry by performing a selection on the proper 
motions between our epoch and the archival ACS GO-11633 
data set (PI: Rich), centered on the same field. 2,430 sources 
were flagged as cluster members, of which 215 belong to the 
white dwarf (WD) sequence (see the CMD on Fig.[Tl left panel). 
The remaining 2,215 MS stars are red dwarfs ranging from MOV 
spectral type {R^ ^ 0.57?©, Mp ^ 0.5Mo) down to the hydrogen- 
burning limi t, as confirmed by superimpo sing on the CMD a n 
isochrone bv lDotter etaH (l2007h . as done in lRicher et an(l2008h . 

3. Systematic correction 

The instrumental magnitude -2.5 log(DN) of each star was reg- 
istered to the median instrumental magnitude of stars measured 
on the deepest frame of the F814W series. We refer to this mag- 



nitude as m. Saturation occurs at m < -13.4. The RMS cr,,, of 
our full sample of 5,078 light curves has been compared with 
the expected noise budget, as calculated by combining theoreti- 
cal photon-, dark-, background- and readout noises (right panel 
of Fig. [1] orange line). Individual cr,„ have been averaged on 0.5 
mag bins with a clipped median (red circles in the same plot) to 
exclude outliers from the comparison. 

On average, the observed noise level on the bright side is 
way higher than the expected, even by 50-60% for stars with 
m < -12 (red circles, right panel of Fig. [T]i. We identified the 
source of most of this excess noise: a variation in the photomet- 
ric zero point (ZP) induced by systematic changes of the PSF 
shape. Long -term instability of the P SF was already reported 
for ACS bv iAnderson & Kind (l2006h . We noticed that this ZP 
change (hereafter called PSF-ZP) follows a well-defined pattern 
as a function of time and average x, y position on the detector, 
measured in physical pixels. The pattern can be mapped by eval- 
uating for each star two diagnostic parameters that appeared to 
be strongly correlated with the PSF-ZP shift: 

1. the difference between the median magnitude (m)beg mea- 
sured on frames taken at the beginning of the orbit, and the 
median magnitude (m)end from frames taken at the end of the 
orbit (left panel of Fig. |2] color-coded in the range -0.02- 
0.02 mag from black to red). 

2. the difference beetwen the median magnitude of the star 
(m)^* measured on 16 consecutive frames taken during the 
seventh "visit" of the program (2453451 < JD < 245352), 
and the median magnitude (m) of its full light curve (right 
panel of Fig. [2] same scale); 
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average x(raw) average x(raw) 

Fig. 2. Mapping the PSF-ZP shift as a function of time and position on the ACS detector, with two different diagnostics (see text). 
First {left panel) as difference between the median magnitudes (m)beg and (m)end measured on frames taken respectively at the 
beginning and at the end of the orbit. Second {right panel) as difference beetwen the median magnitude of the star (m)^* measured 
during the seventh "visit" of the program (2453451 < JD < 245352), and the median magnitude (m) of its full light curve. Top 
panel: All the high-S/N light curves (cr„, < 0.02) registered to their average magnitude (m). In all the panels the color scale spans 
the range -0.02-0.02 mag from black to red). 



The pattern is very similar in both cases. The first diagnostic 
{m)bee - {m)eni is probably a proxy to the real origin of the PSF- 
ZP systematic: a thermal/mechanical instability linked to the or- 
bital phase. The dependence of PSF-ZP on time becomes evi- 
dent when all the light curves with a high S/N (i.e., measured 
on all the 252 frames and having cr,„ < 0.02) are registered to 
their average magnitude (w), stacked in the same plot, and color- 
coded as a function of (m)^* - (m) (Fig.|2] top panel, color scale 
from black to red in the range -0.02 to 0.02 mag). It is clear 
that on average stars whose flux is overestimated during the sev- 
enth visit are also sistematically underestimated in the last visits 
(JD>2453460), and viceversa. 

For a given frame, and within the same chip, the PSF-ZP is a 
smooth function of the position on the detector The diagnostics 
(m)^* - (ot) and (;7j)beg - {m)sni are too noisy when evaluated 
on faint stars to implement an effective correction with them. 
We chose instead to correct the PSF-ZP with a local approach, 
adapted from the differential photometry algorithms described 
in[Nascimbeni et al. (201 1). For each target star / in our sample, 
a set of nearby reference stars k - I, . . .,N has been selected 
with the following criteria: 1) they had to lie on the same chip of 
the target and within 200 pixels from it, 2) they had to be at least 
twenty and their total flux and had to exceed ten times the flux of 
the target, and 3) they had to be detected on at least 250 frames 
over 252, instead of the 200-frames limit required for a target 
star When the requirement 2) was not met, the radius was in- 



creased until needed. If on a given frame j a reference star k was 
not detected, or when its magnitude mji^ was more than 30-^,* 
off its value averaged over the series {nijk), its magnitude was 
set to (mj^k). A reference magnitude niQ j was calculated on each 
frame by performing the weighted mean of the magnitude s of the 
reference stars: moj - Zkimik/crjj)/ Zki^/o-ji) (Broeg et al.1 
2005). The target magnitude was then normalized to niQ j. The 
PSF-ZP correction was applied only when it decreased the over- 
all RMS of the target light curve. After the correction, the me- 
dian RMS of the light curves (right panel of Fig. [1] cyan circles) 
is greatly improved. For the brightest stars {m < -12), it is de- 
creased to a level ~ 15% above the expected noise. 

Most of the remaining exceed noise is due to local ZP 
changes enhanced by the dithering pattern employed: an integer 
Ax, Ay shift plus a small sub-pixel offset 6x, 6y were added to the 
initial pointing XQ,yo (units are in physical pixels). We call this 
systematic effect dith-ZP. Each light curve, already corrected for 
PSF-ZP, underwent two decorrelating algorithms: 

1 . the median magnitude of the star (m)Aj:,Ay was calculated for 
each subset of frames sharing the same integer-pixel dither 
Ajc, Ay. The magnitude m of each frame in the subset Ajc, Ay 
was then registered to (OT)A.v,Ay 

2. for each Ax, Ay subset corrected by the previous step, we 
considered the magnitude nij on each frame j as a function 
of the sub-pixel shift 5xj, Sy j. A fit by ordinary linear least 
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squares was caiTied out to find the coefficients cq, Cv, c,. such 
that 

nij - cq + Cx 6xj + Cy dyj . (1) 

Once the best-fit value m'j = cq 6xj+Cy Syj was evaluated 
for each frame, the corrected light curves were extracted and 
normahzed to zero by evaluating m'j - nij. 

Both steps 1) and 2) were applied only when the RMS de- 
creased. The resulting median RMS of the PSF-ZP + dith-ZP 
correction is plotted on the right panel of Fig. [T]as green circles. 
The small black points on the same plot show the RMS of each 
single light curve of cluster members after our PSF-ZP and dith- 
ZP corrections have been applied. Their noise level approaches 
the theoretical limit on every magnitude bin, demonstrating that 
the algorithms we have applied have been eff'ective. 

4. Light curve analysis 

4.1. Search for transit-like events 

We searched for transits in the full set of 5,078 light curves (that 
is, including fi eld stars) applying t he box-fitting least-square al- 
gorithm (BLS. lKovacs et al.ll2002h . For each curve BLS was set 
to search for periodic dips of duration A and depth 6 with 10,000 
trial periods between 0.2 and 14 days. The relative transit dura- 
tion q - Is. IP was constrained to values possible for planetary 
transits around low-MS stars (7?* = Q.Q^-I.ARq). 

For each star three detection diagnostics have been calcu- 
lated: the signal residue (SR), the signal detection efficiency 
SDE) associated to the maximum peak in the BLS periodogram 
Kovacs et al.ll2002h . and the detection S/N ratio defined as 

S/Nbls = - ■ (2) 
cr 

where cr is the (unbinned) photometric noise and Ux the number 
of data points sampled during transits. 

Other more sophisticated detection diagnostics, such as the 
"signal-to-pink" S/N (Pont et al. 2006; Hartman et al. 2008), are 
robust if correlated noise CTrsd ("red noise' MPont et alJ l2006) IS 
present. However they require the knowledge of cr^^i on time 
scales close to A. This is difficult to be evaluated on our data, 
as transits are expected to be undersampled by the observing 
cadence. Transits of a P ~ 3d planet around a M4V star are 
expected to last 7?*f/(7ra) ~ 60 min at most, and only ~ 30- 
40 min for later types, while images are sampled every 32 and 
64 min. However, the amount of red noise here is very low, as 
demonstrated by the similarity of measured RMS to the theoret- 
ical one. Therefore we decided to employ both S/NgLS SDE 
as detection diagnostics. 

To set a reliable detection criterium, 2,215 light curves were 
simulated, having the same sampling times f, and noise level of 
the real M-dwarfs. A synthetic transit (following the analytical 
model bv lMandel & Agolll2002h of a 1 TJjup planet was injected 
in each curve, with a random uniform distribution in P and sin /. 
P was bounded to the range 1-5 d, while sin / was constrained 
to allow transits. The process was iterated 200 times, for a to- 
tal of 450,000 injections. Then we tried to recover the transits 
with BLS, setting the same parameters used for the real search. 
A planet was defined "recovered" if at least two transits were 
sampled, and if the estimated orbital period (or a low-order har- 
monics: 2:1, 3:1, 3:2) matched the injected one. The distributions 
of the "injected" and "recovered" transits in the SDE vs. S/Nbls 
plane are plotted as black and red points in the upper left panel 



of Fig. [3] By dividing the parameter space in cells and evaluat- 
ing the fraction of "recovered" over "injected" transits (Fig. [3] 
lower left panel), we obtained an estimate of the expected frac- 
tion / of real transits successfully detected ("true positives"). 
We defined detection criteria that guarantee a fraction of false 
positives smaller than 10% (that is, / > 90%): 

( SDE > 5.25 

{ S/Nbls > 65 - 9 ■ SDE (3) 
\ S/Nbls > 7 

(blue line in Fig. [3]). With this choice, the fraction of false nega- 
tives (that is, real planets discarded by selection criteria) is about 
55%, and gets larger for planets smaller than Jupiter This is un- 
avoidable, if one wants to keep as low as possible the fraction of 
false positives. 

The position of all the real sources in the (SDE, SNbls) 
plane is shown in the upper-right panel of Fig. [3] Only four stars 
among the full sample meet the criterium set in Eq. (O or get 
very close to the threshold. These "bordeline" targets (ID#0269, 
1961, 5936 and 7637) were inspected and cross-checked indi- 
vidually. 

ID#269 light curve is crippled by a CCD bad pixel falling 
just under the star in one of the dithering positions. 

ID#1961 is contaminated by brighter sorrounding stars, with 
strong decreases of flux in one fourth of the images. The BLS 
signal is probably spurious. 

ID#5936 seems well measured by the reduction pipeline, 
though it lies extremely close to a saturated star. Its light curve 
should be treated with caution. The parameters of the detected 
signal (P ~ 2.1 d, 6 - 0.08 mag, q = 0.025) would be com- 
patible with a ~ 1 /?jup transiting body with zero impact param- 
eter, or with a grazing eclipsing binary. We classified ID#5936 
as a cluster M-dwarf member analyzing its proper motion (with 
= 0.34Mq, Ri, - 0.32Rq from its color), though its posi- 
tion in the CMD diagram is offset from the MS by 1 .4 mag in 
mF8i4w and 0.2 mag in color (Fig[3] lower right panel). This can 
not be due to binarity alone, and maybe the presence of a bright 
contaminant or the departure from the normal evolution of the 
companion play a role. There is also a non-negligible probability 
that ID#5936 is a field star having a proper motion compatible 
with the common motion of the cluster. It it worth noting that at 
least one data point with a similar decrease in brightness (0.08 
mag, that is ~ 4cr) fell outside the expected transit windows fit- 
ted by BLS, and that the estimated duration (A = 76 min) is way 
larger than the expected. We do not consider ID#5936 a convinc- 
ing planetary candidate. Its coordinates are listed in Table [T] for 
possible further studies. 

ID#7636 was rejected because it fell over a bad column on 
frames corresponding to "transits". 

In summary, we did not detect transits on our light curves, 
at least with an acceptable degree of statistical significance. We 
discuss the significance of this null detection in Section|5] 

4.2. Searcli for variable stars 

We performed a search for variable stars in our full database 
of 5,078 light curves corrected fo r systematic error s. First, the 
coefficient of spectral correlation dFerraz-Melldl 198 1) was cal- 
culated for each Ught curve. Following the method described in 
ide Marchi et a D (120071 l20T(l . we obtained a sample of thirteen 
suspected variable stars (Table [T]i. All these candidates, based 
on their proper motions and position in the CMD, were identi- 
fied as field stars with high confidence. To classify these objects. 
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Fig. 3. Upper left: distribution of SDE and SNbls from the BLS analysis of 450,000 artificially injected transits on synthetic light 
curves. "Recovered" transits are plotted as red dots (see text). Lower left: same as above, where the parameter space has been 
divided in cells and color-coded as a function of the fraction / of transits successfully recovered. Cells with / > 90% (that is, with 
an incidence of false positives smaller than 10%) are highlighted with a black border The blue line correponds to the threshold 
defined in Eq. Upper right: distribution of SDE and SNbls from the BLS analysis of the full sample of 5,078 real light curves. 
Four low-significance candidates are labeled. Lower right: location of the four low-significance candidates on the CMD. 



a leas t-squares iterative sine-wave search was applied (IVanicekl 
Il971h . 

Most of our candidates show a single harmonic sinusoidal 
shape, and short periods (P < 9 d): ID#830, ID#7523, ID#2178, 
ID#5600, ID#4957. Without other elements, it was not pos- 
sible to derive an unambiguous classification for these vari- 
ables. We suspect that these stars are most probably field BY 
Draconis variables, i.e. spotted and rotating KM dwarfs. This 
tentative classification is supported by their very red colors 
('«F606W - '«F814W - 1.27-2.65). 

On four stars a good best-fit can be obtained using two har- 
monics. The second harmonic in the light curves of ID#3428, 
ID#1383, ID#4430 and ID#270 could indicate the presence of 
spots on the surfaces and confirm the BY Dra-type classification, 
while in the light curve ID#6119 the second harmonic reveals 
the profile of a W UMa contact eclipsing binary system. Two 
stars (ID#2086 and ID#258) show evident orbit-to-orbit variabil- 
ity but time coverage is too short to infer reliable values for their 
periods: we classified them as generic "long period variables". 



Finally, the fluctuations in the light curve ID#1882 are too small 
to allow us to confirm its nature of variable star. We discarded 
this last candidate. The summary classification of the entire sam- 
ple is reported in Table [1] along with the best candidate transit 
found by BLS and discussed in Section l4n 

5. Completeness and significance 

The significance of our null detection of transits was assessed by 
considering only the 2,215 cluster- member M dwarfs for which 
we derived reliable estimates of and from their position 
in the CMD. 

For a planet of given radius Rp and orbital period P, the num- 
ber of expected planet detections is given by: 

A^p = JJ (Op(P, Rp) ■ Ogeo ■ (Ddet) dP dRp (4) 

where Op is the fraction of stars with a planet, Ogeo is the ge- 
ometric a priori probability for that system to be aligned such 
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Fig. 4. Top panels: light curves of the variable stars found (first twelve entries in Table[TJ, folded around the best-fit period. Bottom 
left panel: position of the variables (red IDs) in the (mpeoew - 'Mf814w,'Wf8i4w) color-magnitude diagram. Bottom right panel: 
photometric RMS of the variables (red IDs) compared with all the analyzed light curves. 



that a transit occurs, Odet is the probability for that transit to be 
detected by our pipeline, A^* is the number of target stars. As in 
our case we found A^p = 0, we wanted to estimate an upper limit 
to Op, at least in the (P, R^) range for which the efficiency of our 
search Ojet is not negligible. Odet is expected to be dependent 
upon the transit depth (Rp/Ri,)^, the duration A, and the orbital 
period P. 

For simplicity, our analysis was limited to two values of 
planetary radii: "jupiter" planets (1 ^jup) and "neptune" planets 
(0.338 Rjup)- To estimate Ojet, we ran simulations in a way simi- 
lar to what has been done to set the detection threshold (Section 
14.1b . In each set of 2,215 simulated light curves, a synthetic tran- 
sit was injected in each curve, with a random uniform distribu- 
tion in P and sin / (1 < P < 5 d, sin / was constrained to al- 



low transits). The process was iterated 200 times, for a total of 
450,000 injections. Then we tried to recover transits with BLS, 
setting the same parameters and detection criterium defined in 
Eq. (O and adopted for the real search. To derive Odet, the ra- 
tio between detected and injected transits was evaluated for each 
0.1 -day bin of period P. The resulting distribution is plotted as 
a red line in the left panels of Fig. |5] as a function of the in- 
jected orbital period Pin. For comparison, the ratio between "re- 
coverable" transits (i.e., with at least two transits sampled and 
P/Pin =1:1, 3:2, 2:1) and injected transits is plotted with a green 
line on the same panels. As expected, Ojet is a decreasing func- 
tion of P, with minor features at integer and semi-integer values 
of P due to phasing effects. Note that for "neptunes", Ojet is 
extremely low (0.005-0.01) even at short periods (P ~ 1-2 d). 
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Table 1. Parameters for the variable stars found. 



ID# 


RA (2000.0) 


DEC (2000.0) 


"JF814W 




p.m.(Q') 


p.m. ((5) 


P 


Am 


^ harm 


notes 




h:m:s 


d:m:s 






(pixels) 


(pixels) 


(days) 


(mag) 






3428 


17:41:06.078 


-53:45:47.62 


21.884 


1.513 


-1.2192 


0.4788 


2.860 


0.047 


2 


double-wave BY Dra 


1383 


17:41:11.653 


-53:45:02.54 


20.538 


1.848 


-1.3037 


0.2494 


2.650 


0.025 


2 


double-wave BY Dra 


270 


17:41:14.751 


-53:45:04.94 


20.110 


1.213 


-1.0478 


0.7856 


2.290 


0.044 


2 


double-wave BY Dra 


6119 


17:40:58.818 


-53:45:41.10 


20.054 


1.367 


-1.4436 


0.5940 


0.378 


0.051 


2 


WUMa 


4430 


17:41:03.155 


-53:44:49.11 


20.456 


1.761 


-1.2356 


0.4401 


1.095 


0.069 


2 


BY Dra? 


830 


17:41:13.183 


-53:45:37.72 


22.598 


2.654 


-1.3962 


0.6337 


2.700 


0.026 


1 


BY Dra? 


7523 


17:40:55.370 


-53:45:27.03 


20.306 


2.463 


1.0304 


0.5997 


2.520 


0.009 


1 


BY Dra? (weak) 


2178 


17:41:09.636 


-53:43:17.61 


19.362 


2.146 


-1.6238 


0.5057 


2.790 


0.007 


1 


BY Dra? (weak) 


5600 


17:41:00.192 


-53:42:46.28 


19.979 


1.271 


-1.5276 


0.3267 


« 8.9 


0.025 


1 


BY Dra? 


4957 


17:41:01.787 


-53:43:31.37 


20.778 


2.520 


-1.6460 


0.3394 


~ 8 


0.019 


1 


BY Dra? 


2086 


17:41:09.888 


-53:43:10.59 


19.359 


1.097 


-1.5737 


0.3106 


12 


0.006 


1 


long period 


258 


17:41:14.789 


-53:44:58.36 


19.786 


0.813 


-1.7434 


0.3614 


« 25 


0.008 


1 


long period 


1882 


17:41:10.379 


-53:44:52.05 


19.430 


0.792 


-1.0721 


0.7363 


0.920 


0.005 


1 


spurious candidate? 


5936 


17:40:59.316 


-53:44:06.51 


20.524 


1.493 


0.0181 


0.0055 


2.120 


0.080 




grazing binary? 



Notes. The columns give: the ID number of the star, the right ascension a and declination 6 at epoch 2000.0, the calibrated magnitude m(F814W) 
and color ra(F814W)-ra(F606W), the proper motion (a,6) in ACS-WFC pixels relative to the cluster, the period found in days, the amplitude 
found in magnitude, the number of harmonics employed in the fit and a tentative interpretation. 



This is a consequence of inefficient sampHng, which makes the 
significance of neptunian transits very weak: 3-4cr even for the 
most favourable case. 

The geometric factor Ogeo was calculated for each injected 
transit as {Rp + Ri,)/a (a is the semimajor axis), and then con- 
volved with Odet to obtain the probability to detect transits on 
a star which is known to host a planet on a random orbit, as a 
function of its period (Fig.|5] middle panels). 

We para metrized the "p lanet occurence" following the analy- 
sis of Howard et alJ (1201 ih on the distribution of 1,235 planetary 
candidates detected by Kepler. <l>p(f ) has been assumed to be a 
power law modified with an exponential cut-off at period Pcut^ 



dlogf 



(5) 



From [Howard et al.l ( 1201 ll) the following parameters were 
adopted: k = 0.0025, yS = 0.37, Pcm = 1-7 days, y = 4.1 
for "jupiters" (8Re < Rp < 32Rg,), and k = 0.002, /? = 0.79, 
f cut = 2.2 days, y - 4.0 for "neptunes" (4Ris < Rp < SR^,). 
^p{P) is plotted, with an arbitrary normalization, as blue squares 
in the middle panels of Fig. |5] 

We first normalized <l>p(P) imposing '^p(P) — 1' that is 
assuming one planet with 1 d < P < 5 d per star The total 
number of expected detections A^exp(Zi •I'p = 1) within each bin 
over the range 1 < P < 5 d is (Odet ■ '^geo ■ •I'p) Ni, . By summing 
over the range 1 < P < 5 d, we obtained A^exp = 23.8 expected 
detections of "jupiters" and 0.14 of "neptunes" (right panels of 
Fig-El blue symbols). For a flat <l>p(P) distribution, A^exp is larger: 
42.3 and 0.45, respectively. 

The upper limit Op max to ^p(P) suggested by our null de- 
tection can be evaluated by simple binomial statistics, normaliz- 
ing Op in order to get a 68.27% (Icr) or 95.44% (2cr) probabiUty 
of zero detections. We estimated for ju piters On niax(lo") = 4 .8% 
and Op,max(2o-) - 12.9% assuming the lHoward et al.l (1201 lb Op, 
and Op,,nax(lo") = 3.3%, Op,max(2cr) = 9.1% assuming a flat Op. 
As expected. Op ^ax is well above unity for the Neptune sam- 
ple, leaving this planetary population essentially unconstrained 
by our data (Fig.|5] lower right panel). 



6. Discussion and conclusions 

We performed a search for planetary transits and variability 
among 5,078 stars imaged in one of the deepest ACS fields ever 
observed, originally aimed at probing the bottom of the main 
sequence of the metal -poor globular cluster NGC 6397. The 
sample includes 2,215 M0-M9 dwarfs of proven membership. 
Though these data were not optimized for such a study, this is 
the largest homogeneous sample of M dwarfs ever searched for 
variability. 

Instrumental drifts and systematic errors caused by dither- 
ing required a careful empirical correction, described in Section 
|3] We developed and implemented algorithms that allowed us 
to approach the theoretical noise limit on the whole magnitude 
range 19 < mF8i4w ^ 26. The brightest cluster members (MOV) 
were measured with an average scatter of 0.003-0.004 mag over 
a time span of 28 days, demonstrating the power of our decor- 
relating techniques and the feasibility of transit searches in the 
low main sequence of GCs. 

We found no valid planetary transit above the significance 
threshold that we set from simulations. Considering only cluster 
stars, whose physical parameters can be reliabily infered, this 
null detection sets an upper limit to the fraction of stars host- 
ing a P < 5 d Jupiter-sized planet to about Op - 4.8% at 
1-cr confidence and 12.9% at 2-cr, assuming the planetary ra- 
dius distribution derived by iHoward et al.l (1201 ll) from Kepler 
data. In other words, only 0.13 detections are expected assum- 
ing that the underlying planetary population is similar to that 
studied by Kepler Most studies based on RVs also hypotesized 
Op < 1% for short-pe riod, Jupiter-sized planets around solar- 
type stars dMarcv et alj2 005). Furthermore, O p is expected to be 
a very steep function of the stellar metallicity ("Fisch er & Valentil 
2005). Therefore, we are unable to make any firm conclusion 
about the occurrence of giant planets in NGC 6397. 

As demonstrated in Section 14.11 our data set is not enough 
sensitive to Neptune-sized planets to draw any conclusion about 
their occurence, though a much hi gher Op is expecte d for M 
dwarfs by Howard et al. (20lT|) and iLovis et al ] (|2009l). among 
others. This was due to inefficient sampling of the time series 
available from archive material, which translated into poor phase 
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Fig. 5. Completeness tests for cluster stars based on artificially injected transits, for 1 Rj^p planets (upper row) and 0.336 /?jup planets 
(lower row). Left panels: detection efficiency Ojet as function of the input period, for planets potentially recoverable (green symbols) 
and for planets effectively detected by criterium in Eq. (O (red symbols). Mid^/e panels: same as above, but Odet is convolved with 
the geometrical probability Ogeo for a planet to transit. The Howard et all ( |201 1|) <l>p(P) period distribution function is plotted in blue 
symbols (arbitrary normalization). Right panels: number of expected transit detec tions per period bin , assuming one planet per star 
within 1 < f < 5 d. 3>p(f ) is assumed to be flat (red symbols) or as modeled bv iHoward et al.l (1201 1.) (blue symbols). 1- and 2-cr 
upper limits for the planet occurence <l>p,max are shown. 



coverage and severe undersampling of transit-like events, whose 
duration is expected to be of the same order of the effective ca- 
dence. 

Twelve new variable stars have been identified in the 
NGC 6397 field (Table E). Most of them can be classi- 
fied as BY Draconis variables, that is, spotted rotating KM 
dwarfs. Interestingly, no variable has been detected among the 
2,430 cluster members. Hundreds of member early-M dwarfs 
(^MpSMw < 21, Fig.[T]) follow the expected noise on time scales 
up to 28 days, though they were measured with a 0.003-0.006 
mag precision. The lack of eclipsing binaries is not surprising. 
The number of expected detections car i be estimated by scaling 
down the number of EBs detected bv lAlbrow et alJ (1200 1) for 
47 Tuc. If one takes into account the lower fraction of binaries 
in NGC 6397 (< 3%, even at radii smaller than the half-mass 
radius, Milone et al ] 120111) and the smaller number of targets 
(2,215 vs. 46,422), we would expect much less than one EE. On 
the other hand, the lack of BY Dra variables is more puzzling, as 
one would expect 4-5 such detections, considering the number 
of monitored clust er stars in this paper . We note that our stars 
are cooler than the lAlbrow et al.l (1200 ll) sample, and we should 



expect longer photometric periods. Our search is not sensitive to 
periods > 28 days, so this could be a possible explanation. 

We evaluated an upper Umit to the photometric jitter o-jit of 
the brightest members (MOV) by subtracting the contribution of 
the expected noise os^p to the measured scatter o-gbs (that is, as- 
suming cr?j - cr^^^ - cTgjjp). The fraction of stars / with o-jit > 2 
mmag is / < 2%. This va l ue sho uld be compared with the re- 
sults found bv lCiardi et alJ (1201 ll) examining the first quarter of 
Kepler photometry on 2,182 field M dwarfs: these data cover 
an interval of 33 days with an average cadence of 30 minutes, 
quite similar to our cadence and time scale. They found a frac- 
tion / ^ 20% of stars with cr > 2 mmag, that is at least an order 
of magnitude larger fraction than what we measured in NGC 
6397. The low MS of this cluster is extremely stable and there- 
fore worth targeting using more optimized observations. JWST 
NIRCam, for instance, would be able to probe the bottom of the 
MS of NGC 6397 (mpsM-w ~ 24) with a photometric precision 
better than 0.01 mag on a single 600 s exposure, without any of 
the coverage/sampling issues mentioned above. 

Most of the analysis techniques presented in this Paper can 
also be applied (with little or no modification) to other exist- 
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ing ACS/WFC3 time series of rich stellar fields. This is the case 
for the metal-rich globular cluster 47 Tucanae, that has been 
imaged with ACS and WFC3 with a larger coverage. Such a 
search for tran sits will allow us to complement the results by 
iGilUland et all (I2OOQ) in a diff'erent range of spectral types and 
planetary masses. 
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